Scale dependence of twist-three contributions to single spin asymmetries 
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We reexamine the scale dependence of twist-three correlation functions relevant for the single 
transverse spin asymmetry in the framework of collinear factorization. Evolution equations are 
derived for both the flavor-nonsinglet and flavor-singlet distributions and arbitrary parton momenta. 
Our results do not agree with the recent calculations of the evolution in the limit of vanishing gluon 
momentum. Possible sources for this discrepancy are identified. 
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I. INTRODUCTION 

Large transverse single spin asymmetries (SSAs) have 
been observed in different hadronic reactions and these 
observations generated a lot of interest. Such exper- 
iments are conceptually rather simple, but their the- 
oretical description proved to be challenging as the 
leading-twist contributions to such asymmetries vanish, 
see [Hjlil El for a review. Over the past few years there 
was a splash of theoretical activity in this field, which 
mainly followed two lines: the /cj_ factorization in terms 
of the transverse-momentum dependent (TMD) distri- 
butions (e.g. |, 0, H, 0, @, S M E3). or > alterna- 
tively, collinear factorization including twist-three con- 
tributions in terms of multiparton correlation functions 
[H 0, [3 [H, M EE El, 111 ID, [U, Hp. These two 
techniques have their own domain of validity and were 
shown to be consistent with each other in the kinematic 
regime where they both apply [23|, [HI, HE HE] . 

However, practically all existing calculations have been 
so far at the leading order (LO) which corresponds, 
roughly speaking, to the (generalized) parton model. In 
order to test the QCD dynamics and eventually also re- 
duce the dependence of theory predictions on both the 
factorization scale and the renormalization scale of the 
strong coupling, it is necessary to calculate the scale 
dependence (evolution) of the relevant nonperturbative 
functions. The corresponding calculation was done in 
Refs. JU [HI where the evolution equation was derived 
for the gluonic pole contributions to twist-three correla- 
tion functions that are relevant to SSA. Unfortunately, 
it turns out that the equations derived in [U [H[ are in 
contradiction to earlier results obtained in the framework 
of the operator product expansion (OPE). A resolution 
of this discrepancy presents the main motivation for this 
study. 

Although we do reproduce the terms that appear in 
Ref. [27], we obtain two additional contributions to the 
evolution equation both for flavor-nonsinget and flavor- 
singlet operators. We have found that part of the dis- 
agreement is due to omission of cut diagrams correspond- 
ing to the mixing of conventional twist-three correla- 



tion functions that are discussed in the context of SSA 
with the contributions involving a gluon and the quark- 
antiquark pair on the opposite sides of the cut. These 
are the same contributions that have been discussed for 
some time for semi-inclusive DIS [H [2(| and the 
new observation in this work is that they necessarily con- 
tribute to the leading logarithmic accuracy through the 
(factorization) scale dependence. For another part we do 
not have an intuitive explanation. Our OPE-based cal- 
culation contains additional contributions (compared to 
(27l |28| ) for which sending the gluon momentum to zero 
has to be done with caution, as will be explained in the 
text. 

Another question that we want to address is to clar- 
ify the relation of these results to earlier calculations of 
the scale dependence of twist-three correlation functions 
that contribute to inclusive reactions, pioneered by the 
study of the structure function g 2 (x,Q 2 ) in Ref. [3C| . 
We would like to emphasize that the scale dependence 
of an arbitrary twist-three operator in QCD, hence ar- 
bitrary twist-three light-cone correlation function can be 
determined in terms of the two-particle evolution kernels 
introduced by Bukhvostov, Frolov, Lipatov and Kuraev 
(BFLK) [3l|. The BFLK approach has become stan- 
dard in calculations of the spectrum of anomalous di- 
mensions (dilatation operator) in supersymmetric theo- 
ries that are relevant to the AdS/CFT correspondence 
(see e.g. [HI, HI]), but, unfortunately, remains to be 
largely unknown to the broad QCD community. We will 
show that evolution equations for the particular parton 
distributions relevant for SSA (27j do not require an inde- 
pendent calculation but can be obtained from the BFLK 
kernels by simple algebra. Using an updated version of 
this technique [34l . l3q | we derive the complete evolution 
equations for all relevant three-particle correlation func- 
tions for arbitrary gluon momenta. This generalization 
is necessary since the gluon-pole contributions considered 
in [27], HI] do not have autonomous scale dependence. 

The presentation is organized as follows. Section 2 is 
introductory. We review basic properties of three-particle 
correlation functions and introduce a very convenient 
decomposition of momenta (coordinates) and the field 
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T^(zi,Z2, z 3 ) = gq(z 1 n)-f + F l _ l+ (z 2 n)q(z 3 n) , 
AT ll (zi,Z2, z 3 ) = gq(z 1 n)'y + ^ 5 iF^ + (z2n)q(z 3 n) (1) 



operators in the spinor representation. This rewriting trix elements of the quark-antiquark-gluon light-ray op- 
makes the symmetries explicit and drastically simplifies erators 
the forthcoming algebra. Section 3 contains a detailed 
derivation of the flavor-nonsinglet evolution equation in 
the BFLK approach and the comparison with [2?], [28|. 
Our result coincides identically with the corresponding 
equation for the quark-antiquark-gluon correlation func- 
tion relevant for the structure function q 2 (x,Q 2 ). The 

&*>\ T ~* / — ( f d) ~~ ( f d) 

flavor-singlet evolution is considered in Section 4. In this and also the correlation functions Tg F and cor- 
case there are two independent evolution equations for responding to the three-gluon operators 
positive and negative C-parity that involve three-gluon 

correlation functions involving SU{2>) structure constants , syabcT?a / \j?b i \t?c i 

f abc and d abc , respectively. The equation for positive C- G n P \\ z ^ z ^ z ^> = # 6 ± F + P ^ n ) F +^ n ) F +\\ z ^ n ) ■ 
parity coincides, again, with the corresponding equation (2) 
for the structure function g 2 (x,Q 2 ), whereas the equa- 
tion for negative C-parity is, to our knowledge, a new 
result (in this form). Finally, in Section 5 we summa- 
rize. The full list of the BFLK kernels in the momentum 
representation is collected in the Appendix. 



Here is the light-like vector, n 2 = 0, the "plus" pro- 
jection is defined as a + — a M n M . Note that there are 
two gluon operators with a different color structure; the 
factors C± c are written in terms of the SU(3) structure 
constants 



II. GENERAL DISCUSSION 

A. Definition and support properties of 
three-particle light-cone correlation functions 



(3) 



Following Ref. [27[ we will consider the correlation 
functions T q ^ F , Taq,f, corresponding to the nucleon ma- 

I 



In all cases the path-ordered Wilson lines are implied that 
ensure gauge invariance. They are not shown for brevity. 



The definitions of parton distributions used by Kang and Qiu are (cf. Eqs. (12), (16), (14), (23) in [271]): 

1 r dz 1 dz 2 jp + ( Zl x+z 2 x 2 



T q . F {x,x + x 2 ) 



f^f(x,x + x 2 ) 



(2tt) 2 
1 f dz?dz. 



(P,aT|a M T M (0,«2,«i)|P,« r ), 



/ ^0e iP ^*+*^ (P, s T \s»AT»(0, z 2 , Z1 )\P, st) , 



(27T) 



Ti^ F (x,x + x 2 ) = *j0 e ir + fr*+»*> )(P,s T \s»G% x (0 ) z 2> z 1 )\P,s T ). (4) 



Here s p is the nucleon spin vector normalized by the condition s 2 = — 1, s M = —^ ivpa s v n p h a with h being the second 
light-like vector, fi 2 = 0, nn — 1 and = — e P An ™ j 60123 = 1- The spin vector s M is assumed to be transverse, 
ns — ns — 0. In the last two equations T^' and correspond to matrix elements of G+ pA and G~ x , respectively, 
cf. Eq. ©. 

For our purposes it is convenient to use a more symmetric notation with quark, antiquark and gluon momentum 
fractions treated equally 

(P, St^T^zx, Z2, z 3 )\P, s T ) = 2Pl J Vx e - lP +^» x ^ T qFq {x u x 2 , x 3 ) , 
(P, axis' 1 AT^zu Z2, z 3 )\P, s T ) = 2P^ J Vx e - lP +^ AT qFq {x u x 2 ,x 3 ) , 



gP x {P, ST \^G% x {z ll z 2 ,z 3 )\P, ST ) = P% J Vxe- lP +^ x ^T± F {x llX2 ,x 3 ), 
e p ± x {P,sx\s^ pX (z 1 ,z 2 ,z 3 )\P,sx)=Pl J Vxe~ iP +^^ ATi F (x 1 ,x 2 ,x 3 ), 



(5) 
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J T>x — J dxidx 2 dx3 S(xi + x 2 + x 3 ) . (6) 

We assume here the standard (relativistic) normaliza- 
tion of states; the factor P£ is necessary to ensure the 
reparametrization invariance to the choice of the light- 
cone vector n„ — > an^. Obviously 

T q ^ F (x,x + x 2 ) = T qFq (-x - x 2 ,x 2 ,x) , (7) 

and similarly for the other distributions. Writing the 
definition in such a form (Fourier-transformed) makes ex- 
plicit the support properties of the correlation functions: 
they are defined on the surface x\ + x 2 + X3 — and 
effectively are functions of two variables only. 

The support of parton correlation functions in the no- 
tation by Kang and Qiu, i.e. in the (a;, x 2 ) plane, is shown 
in Fig. Q] It can be separated in six different regions 
where x, x 2 and x -\- x 2 are positive or negative, respec- 
tively. The parton-model interpretation of each region 
is different. As discussed in detail in Ref. [36[, interpre- 
tation of light-cone correlation functions as describing 
sequential emission or absorption of quark-, antiquark- 
and gluon-partons by the target arises by choice of a 
particular representation in terms of the sum over in- 
termediate states (cut diagrams) that does not involve 
semidisconnected contributions. In particular, the upper- 
right region in Fig. [T] corresponds to emission of a pair 
of a quark-parton and a gluon with momentum fractions 
x > 0, x 2 > 0, respectively, and subsequent absorption 
of the quark-parton with x + x 2 . 

The picture can be made more symmetric going over 
to the correlation functions (fT]) and using analogue to 
barycentric coordinates (3?| as shown in Fig. [5J 

x = x\&\ + x 2 e 2 + X3S3 = X1E1 + x 2 E 2 . 

The six different regions, labeled (12) + 3", 2+(13) - , 
etc., correspond to different subprocesses at the parton 
level [36j]; For each parton k = 1,2,3 "plus" stands for 
emission (xj. > 0) and "minus" for absorption (xk < 0). 



Alternatively, one may think of "plus" and "minus" la- 
bels as indicating whether the corresponding parton ap- 
pears in the direct or the final amplitude in the cut dia- 
gram 38]. 

It is important that different regions do not have au- 
tonomous scale dependence; they "talk" to each other 
and get mixed under the evolution. The particular mix- 
ing pattern can be understood as follows. 

It is easy to see that to one loop accuracy any con- 
tribution to the evolution equation only involves two of 
the three partons. For example, parton-1 and parton-2 
can exchange a gluon whereas parton-3 stays a spectator, 
etc. The full evolution kernel can, therefore, be split in 
three two-particle kernels that involve parton pairs (12), 
(23) and (31), respectively. Schematically 

H-T qFq = ^[H 12 + H 23 + H 31 ] <g> T qFq . (8) 
[i In 

Each two-particle kernel is, obviously, a function of the 
contributing parton momentum fractions, e.g. H\ 2 = 
) and due to energy conservation x\ + 
x' 2 = x\ + x 2 . In other words, the rate for the scale 
variation of a three-particle parton correlation function 
with given values of momentum fractions x±, x 2 can only 
involve this function on the line of constant x\ +x 2 = —X3 
(for the contribution of Hi 2 ), as shown in Fig. [2] (thick 
blue line AB). 

Since 13^0 (in general), the total momentum frac- 
tion carried by the two participating partons is non-zero, 
xi + x 2 0. This situation is familiar from studies of 
the scale dependence of leadin g-tw ist generalized parton 
distributions (GPDs) (see e.g. [39l. l4fj| ) and the three re- 
gions 2+ (13)-, (12)+3-, l+(23)- traversed by the thick 
blue line AB in Fig. [2] are in one-to-one correspondence 
to the two DGLAP regions (2+1" and 1+2") and the 
ERBL (central) region (l + 2 + ) in the corresponding evo- 
lution equations. As well known [39|, the scale de- 
pendence of GPDs in the DGLAP regions is autonomous, 
whereas in the ERBL mode there are also terms describ- 
ing "leakage" from the DGLAP regions. In the present 
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context, this result implies that evolution equation for a 
generic three-particle light-cone correlation function for 
the momentum fractions X\,X2 in the (12) + 3~ region, as 
in Fig. [21 will receive nontrivial contributions from the 
2+(13)-, l+(23)-, 2- (13)+ and 1"(23)+ regions as well, 
which have a different partonic interpretation. 

So far we have considered the contribution of H±2 only. 
The contributions of H23 and H\ 3 are, in turn, kinemat- 
ically constrained to the lines of constant and 
x 'i + x 3i respectively, and, for the particular choice of 
x\,x% in the (12) + 3~ region, correspond to the DGLAP 
modes of the corresponding GPD-like evolution equa- 
tions. Recall that the DGLAP evolution is ordered in 
the momentum fraction; hence only parts of the kine- 
matically allowed regions contribute, as shown in Fig. [2] 
by the thick red CD (H 13 ) and thick green EF (H 23 ) 
lines. 

Note that evolution of the parton correlation function 
for given values of x%, x 2 only involves the regions of mo- 
mentum fractions outwards from the center. If one draws 
a small(er) hexagon (see Fig. [2]) inside the big one, then 
evolution of the correlation function outside of the small 
hexagon does not depend on the latter inside the small 
hexagon. This "radial" ordering is a generalization of the 
usual momentum fraction ordering in DGLAP equations. 
In contrast, there is no ordering/restrictions in the "az- 
imuthal" direction and, in principle, all regions "talk" to 
each other. This property may result in the increase of 
the evolution rate in the small x region. 



B. Spinor representation and symmetry properties 

Analysis of symmetry properties and the scale depen- 
dence of the light-cone correlation functions can be sim- 
plified significantly by going over to the spinor represen- 
tation. Each covariant four-vector x^ is mapped to a 
hcrmitian 2x2 matrix x: 



Finally, the gluon strength tensor i*) t „ and its dual 
can be decomposed as 



F 



a0,a0 



iF a f3,a/3 = 2 ( e apf<*P + ^/d/j) 



(11) 



Here f a p and f & g are chiral and antichiral symmetric 
tensors, /* = /, which belong to (1,0) and (0, 1) repre- 
sentations of the Lorenz group, respectively. The anti- 
symmetric tensors, e a p — e a and e & g = e a/3 are normal- 
ized by £12 = — £j2 = 1 an d use d f° r rising and lowering 
indices 



u a = e al3 u . 



Ua 
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C a/3 U ' ■ (I 2 ) 

The contractions (ab) and (ab) are defined as 

(ab) = a a b a = -a a b a , (06) = ajf 1 = -a & h- (13) 



The scalar product of two four- vectors a and b takes the 
form a^¥ l = \a ai yb aa . 

For convenience, we present the expressions for Dirac 
matrices in the spinor basis: 



r 







J a/3 







75 



< 
6f 



(14) 



where 75 = 17 7 7 7 . More expressions and various 
useful identities can be found in (34 . l35j . 

A covariant generalization of the decomposition in 
terms of light-cone coordinates and a transverse plane 
(|9|) can be found by observing that any light-like vector 
can be represented as a product of two spinors. In par- 
ticular one can parameterize the light-like vectors n and 
h as follows 



^ad — A a A,^ 



(15) 



x aa =x^) aa , 



where — (1, a), = (1, — <?) and a are the usual 
Pauli matrices. In components 



where A = A^, /x = iv . The standard convention (nn) = 1 
corresponds to the normalization of auxiliary spinors 



(fxX) = (A/2) = y/2. 



(16) 



x Q + x a X\ — ix 2 
X\ + 1x2 x — x 3 



(9) 



so that instead of the usual components x^ = {xq, —x} 
each four- vector is described by its light-cone coordinates 
x± = xq ± x 3 and two complex coordinates in the trans- 
verse plane x\ ± 1x2- 

The Dirac (quark) spinor q is written as 



' 



(X^d): 



(10) 



where rp a , are two-component Weyl spinors, ip a 



The basis vectors in the plane transverse to n, n can be 
chosen X a and XaHa and an arbitrary four-vector x 
represented as 



z Aq-Aq, -j- z fi a 

fJ,a + W A a /id + 10 HaX 



(17) 



where z and z are real and w, w = w* complex coordi- 
nates in the two light-like directions and the transverse 
plane, respectively, cf. Eq. Q. In particular the spin 



vectors s M and s M take the form 



--|A Q /XaS M x - MaAdSA/i| , (18) 
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where s„^ = fi a s a a^ a "= (/UsA). 

The "+" fields (" good components" ) are defined as the 
projections onto A: 



i> + = \ a l(j a , x+ = A Q x Q , /++ 

■0+ = \ a 4>a , X+ = A"XA , /+ + 



The scale dependence of the correlation functions ([4|, 
(|5"|) is determined by the renormalization properties of 
the corresponding light-ray operators |T]), ((21). In order 
to make a connection to the existing results, it is conve- 
nient to go over, for quark-antiquark-gluon operators, to 
another basis [4lf 



Sf(z) = gq( Zl )[F p+ (z 2 ) ±i l5 F p+ (z 2 )] 1+ q(z 3 ) . (19) 

Here and below we use a shorthand notation q(z 3 ) = 
q(z 3 n) etc., and also Sp^(z) = Sp^{z\, z 2 , z 3 ). The opera- 
tors contribute to the structure function g2(x, Q 2 ) for 
polarized deep-inelastic scattering and their renormaliza- 
tion properties are studied in much detail jU, 0, [3, 53, 
EE HE S3 - It is easy to check that 



sPT p (z) = -s" {S+(z) + S7(z)) 



s p AT p (z) 



s" (S+(z)-S7(z)) 



(20) 



The second identity is a consequence of the known repre- 
sentation for the product of two antisymmetric e-tensors 
as the determinant of the 4x4 matrix in terms of the 
metric tensors g^ u . 

In the spinor basis one finds for the operators = 
^t a aa tne following expressions: 

S ta( Z ) = 2 9 ^a^+(zi)f +a (z 2 )tp + {z 3 ) 

+ KX+{zi)f+a{z2)X+{z?,) 
S aa( Z ) = 2 9 ^aX+(zi)f+a(z2)x+(z 3 ) 

+\ a # + (z 1 )f + *(z 2 )il> + (z a j\. (21) 

Taking into account transformation properties of the 
quark and gluon fields under charge conjugation, 



CipC^ 1 = x , 

cfc- 1 = -f T , 



cfc- 1 = -f T , 



where / = f a t a , one derives 

CS*{zi, z 2 , z 3 )C 1 = S p (z 3 , z 2 , z{). 



(22) 



(23) 



Since operators of different C-parity do not mix under 
renormalization, it is convenient to introduce the C-even 
and C-odd combinations 

e ± {z) = S+(zi,Z2,z 3 )±S-(z 3 ,z 2 ,zi), (24) 



where S 1 * = ^S 1 ^. Note that these are not the same 
combinations as the ones appearing in Eq. (|20|) since the 
quark and antiquark field coordinates are interchanged. 
The expressions for the operators 6 ± in terms of chiral 
fields read 

e ± (z) = -i| {s^iz) - sxp&iz)} , (25) 
where 

Q ± (z) = $ + (z 1 )f ++ (z2)il) + (z 3 ) ± X+{z?,)f++{z2)x+{zi) , 
Q ± { z ) = X+(zi)f++(z2)x+(z 3 ) ± + (z 3 )/ ++ (z 2 )V'+(2i) • 

It is easy to see that Q ± (z) = ±[Q ± (z)] t , so that the C- 
even ( "plus" ) and C-odd ( "minus" ) ©-operators are her- 
mitian and antihermitian, respectively, i.e. (& ziz (z))^ = 

±e ± ( z ). 

The operators 6^ have autonomous evolution for the 
flavor-nonsinglet sector, whereas in the singlet sector 
they mix with three-gluon operators of the same C- 
parity: 

j*{z) = -^cf c { s ^n + (z 1 )f ++ (z 2 )r ++ (z 3 ) 

-*Ap/? + (*i)/t + (^)/^+(2a)}.(26) 

The color factors C± c are given in Eq. . In the vector 
notation this definition corresponds to 

J^(z) = 2gCf c S p (1 T P23 ± P12) 

xF + ^( Zl )F b +p (z2)FU z 3), (27) 

where P 23 and P12 are the permutation operators acting 
on the field coordinates, e.g. 

P 23 F + v > a (z 1 )Fl p (z 2 )F* u (z 3 ) = F^ a {z 1 )F b +p {z 3 )Fl v {z2) 

(28) 

etc. 

Taking the nucleon matrix elements of the operators 
e ± and T± one obtains the corresponding parton corre- 
lation functions in momentum space 

(P.arl^OOI-P.'r) = 2P| J Vxe- iP +^ XkZk J^{x). 

(29) 

Here x = {x\, x 2 , x 3 }, the integration measure T>x is de- 
fined in Eq. ©. We use the same notation for the correla- 
tion functions as for the corresponding operators which, 
hopefully, will not result in a confusion. Note that al- 
though definitions of the operators &^ and ^ r± , Eqs. ((25]) 
and (|26[) . involve the spin- vector s p , the dependence on 
s p actually drops out in the matrix elements (|2"9l . 
Hermiticity and C-invariance imply that (cf. [17| ) 

(e ± (x)y = ±e ± (~x), e ± {x) = ±e ± {-x), 

(^ ± (.t))*= T ^ ± (-x), F ± {x)=tF ± {-x). (30) 
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This means, in particular, that the correlation functions 
S ± (x) and J r± (x) are real functions. Notice also that the 
function T~ is symmetric and T + antisymmetric under 
the interchange of the last two arguments: 

J r± (x 1 ,x 2 ,x 3 ) = ^J r± (x i ,x 3 ,x 2 ) . 

The six correlation functions defined in Eq. ([5]) (or, equiv- 
alently, Eq. (H))) can be written in terms of two inde- 
pendent quark-antiquark-gluon functions &^(x) and two 
three-gluon functions T (x) as follows: 



TgFgix) = - [(1 + Pl 3 )& + (x) + (1 - Pl 3 )e~ (x) 

AT qFq (x) = -~ [(1 - Pi 3 )6+(x) + (1 + Pi S )6-(x) 



AT 3 ± F (x)=-±(l±P l3 )f ± (x). 



(31) 



Here T qFq (x) = T qFq (xx,x 2 ,x 3 ), etc., and Pik are the 
permutation operators for the corresponding momentum 
fractions, e.g. P\ 2 T ± (x\, X2, X3) = J 7 (x 2 ,xi,x 3 ). As 
follows from ([3Tj) the functions T^p(x) and AT^ F (x) are 
not independent, 

AJ^^!,^,^) =± T^ F (xi, x 3 , x 2 ) - T 3 ± F (x 2 , xi,x z ) 

(32) 

It follows from ([3H)l and (|3"Tj) that the correlation func- 
tions satisfy the following symmetry relation 



where H is an integral operator. The explicit expres- 
sion for H (at one-loop) can be restored from the cor- 
responding result for S ± in (Hj (see also pi I4HL l47j|). 
Expanding this equation at short distances Zi — > one 
reproduces the mixing matrix for the twist-three local 
operators \M, 0, 51, El • 

Alternatively, the answer for HI can be obtained by sim- 
ple algebra from the known expressions for the one-loop 
two-particle kernels [3l|, [HI]. This technique is general 
and applicable to arbitrary twist-three (and twist-four 
[35I ]) evolution equations. In the next Section we will 
use this approach for a more complicated case of flavor- 
singlet operators so it makes sense to explain the details 
on the present (simpler) example. 

The starting observation is that to one-loop accuracy 
any contribution to the evolution can only involve two 
partons; hence H can be represented as a sum of two- 
particle kernels. Schematically 



h = t?k{ WmA + t+fl+^l + WftM } 

(36) 

where the contractions correspond to the sum of rele- 
vant Feynman diagrams (in light-cone gauge) . The corre- 
sponding expressions were derived originally by Bukhvos- 
tov, Frolov, Lipatov and Kuraev (BFLK) [3l|. The com- 
plete list of the BFLK kernels for arbitrary chiral fields 
is given in Ref. [35[ so that it only remains to contract 
the color indices. For the reader's convenience we collect 
all the kernels in Table U in the Appendix. 
After a simple algebra one gets 



TqFq(xi,X 2 ,X 3 ) = Tq F q(-X 3 ,-X 2 ,-Xi) , 
ATg Fq (Xl,X 2 ,X 3 ) = -ATq Fq (-X 3 ,-X 2 ,-Xl) 

T 3F (xi,x 2 ,x 3 ) = T 3F (-x 3 ,-x 2 ,-xi) , 
AT 3F (xi,x 2 ,x 3 ) = -AT 3 ± F (-x 3 ,-x 2l -xi) . 



(33) 



Authors of [27j | also introduce symmetrized quark- 
antiquark-gluon parton distributions 



T q . F (x,x') = i (T q . F (x : x') + T q . F (x',x)^j , 



Aq,F 



(X, X') = - (TAq,F(x, X') ~ T A q tF (x\ x)J . 



(34) 



As follows from Eq. (|33|) such a symmetrization is 
not necessary since T q ^ F (x,x') = T q ^ F (x',x) and 
AT QtF (x, x') = — AT qtF (x', x). In the expressions given 
below we drop the "tilde" notation for these functions. 



III. FLAVOR-NONSINGLET EVOLUTION 

The flavor-nonsinglet light-ray operators 6 + and 6~ 
satisfy the same evolution equation, 



d 



d 



e ± = 0. 



(35) 



N C W - — Hi-3Cjp, 



(37) 



where, in notation of Ref. |3E 



W12 

W13 



H23 
rl 13 



2H+ , 

P 23 H e 2 j 1] + 2H U . 



(38) 



Here Tiik are two-particle integral operators that act on 
the light-cone coordinates of the i-ih and fc-th partons: 

\Hik<i>]{zi,Zk) = / dz'idz'kHizi, z k \z[, z' k ) 4>{z[, z' k ) . (39) 



These kernels are SX(2)-invariant and depend on the con- 
formal spins of partons that they are acting on. The 
corresponding values are j ' = 1 for quarks and j = 3/2 
for gluons (for the "plus" components) so one has to use 
ji = 1, j 2 = 3/2, j 3 = 1. One obtains, for example, 



/> 1 />a 

[Ht 3 & + ](zx,z 2) z 3 ) = da d(3&- 
Jo Jo 



(z" 3 , z 2 , z 31 ) , 



where 



ik = ZiOt + z k a . 



a = 1 — a . 
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Taking the nucleon matrix element of Eq. ([35)1 one ob- 
tains the evolution equation for the corresponding par- 
ton distribution function. Technically, this corresponds 
to going over from coordinate to the momentum frac- 
tion space {zi, Z2, Z3} — * {x%, x%, X3}. Thanks to energy 
conservation, two-particle kernels in momentum fraction 
space can be written in the following generic form 



where Vx' = dx' i dx' k S(xi + Xk — x\ — x' k ). It is assumed 
that restrictions on integration regions over x\ , x' 2 come 
from support properties of the kernels and the parton 
distributions; the corresponding expressions are collected 
in the Appendix. 



[Hik<p\(xi,x k ) = / Vx' Uixi.x^x'^x'^ipix'^x'k) . 



(40) 



The last step, using the first two relations in Eq. 
we obtain 



\i-T qFq (x) 
A* 

(J,-AT qFq (x) 



4tt V 



P 13 MP 13 ^T qFq (x) + ^-(l 
AT qFq {x) + ^ 



-P 13 UP 13 )AT qFq (x) 

i-p 13 ip 13 )v,(i) 



(41) 



r 




FIG. 3: Integration regions in Eq. (|43 



which is our final result. The "Hamiltonian" EI is de- 



fined in Eqs. (|37| , (|38|) . Explicit expressions for the two- 
particle kernels 7^, Wi, and Ti.~ k are given in the 
Appendix. Note that when applying two-particle kernels 
to a three-particle parton distribution one has to treat 
the latter as a function of three independent variables; 
the condition x± + 22 + £3 = is applied afterwards. 

To compare our result with the calculation in Ref. [27j 
we have to change notation to 

T q , F (x,x') = T qFq (-x',x' -x,x), 
T Aq , F (x, x') ee AT qFq (-x', x 1 - x, x) , (42) 

which are (anti)symmetric functions under permutation 
x «-> x' , T q}F (x,x') = Tq jF (x',x) and AT q , F (x, x') = 
—AT q ^ F {x\x). Using explicit expressions for the two- 
particle kernels and taking the gluon-pole limit x' = x 
(x2 — 0) we obtain (it is assumed that x > 0) 



^Tq,F{x,x) = — \ I — 



-T I ./ , £ 



p f yr ft*, N c f (l + z)T^ F (x,0-(l + z 2 )T q ^ F (^Q 
P gq (z)T q>F (Z,£) + — I H - H - TA q ,F(x,0 



N c T q , F (x, x) 



1 

2N, 



(1 - 2z)T q . F {x,x - £) - TAq,F(x,X- £) 

I 



(43) 



where z = x/£, 



P qq (z) = C F 



1 + z ' 2 3 W1 



and 



(1-*)h 



1-z 



(44) The partonic interpretation of different contributions 

in Eq. (|4"3")l is illustrated in Fig. [3] Note that the condi- 
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2 1 3 

FIG. 4: The exchange diagram 



tion of zero gluon momentum X2 — corresponds to the 
choice of x on the (positive) horizontal axis. The first 
term ~ 7^f(£,£) on the r.h.s. of Eq. P^l corresponds 
to the integration over the horizontal line segment xD 
shown in red in Fig. [3l The terms in T qtF (x,^) are due 
to the integration along the blue, AB, and green, EF, 
lines in the regions l + (23)~ and (12) + 3~, and the ones 
in T qjF (x,x — £) correspond to the contributions along 
the same lines in the regions 2 + (13)~ and (13) + 2 _ , re- 
spectively. By construction, the function T qtF (x,x') is 
symmetric and 7a,j,f(x,i:') antisymmetric under the in- 
terchange of arguments, x <-» x' . In Fig. [3] this corre- 
sponds to a reflection symmetry around the horizontal 
axis. 

The evolution equation (|43|) differs from the corre- 
sponding result by Kang and Qiu (see Eq. (99) in (27| ) 
by the two extra terms in the second line. The last 
term, proportional to l/N c , originates from the kernels 

P23H23 and H12 in Eq. (f38|) . On the diagrammatic 
level this contribution corresponds to the "exchange" di- 
agrams of the type shown in Fig. 3] which correspond to 
mixing of the regions that have different partonic inter- 
pretation: 2+ (13)- l + (23)", or (13)+2 _ <-> (12)+3~. 
Similar contributions have been discussed in a somewhat 
different context in Ref. [29l |. 

The other difference is the extra term —N c T q ^ F (x,x) 
in the second line in Eq. (|43[) for which we do not see any 
obvious explanation. The origin of this term in our cal- 
culation can be traced to the kernels TL12 and 7^23- E.g. 
the second term in the expression for 7ii2, Eq. (|A.18[) . 
gives 



d 

^d/x X2/ 



■in 



X2<lx' 2 



x' 2 (x' 2 - X 2 ) 

X2 



T q>F (x,x + x 2 ) r Tq iF (x, x + x' 2 ) . (45) 



Because of an overall X2 factor, it is tempting to put this 
contribution to zero in the X2 — ► limit. However, it is 
easy to see that in the same limit the integral becomes 
linearly divergent so that at the end a finite contribution 
arises. This term can easily be missed if the gluon mo- 
mentum is put to zero at the beginning of the calculation. 



It is instructive to analyze Eq. (143]) in the large-iV c 
limit. Neglecting 1/N C terms in (|43l) one is left with 



{Pqq(z) + z)T q , F (x,£) -TAq,F(x,0 



(46) 



where P qq (z) is obtained from P qq (z), Eq. ([44]) . omit- 
ting the C F factor. One sees that the scale dependence 
of T q<F (x,x) is determined by T q ^ F (x,^) in the region 
x < £ < 1, which corresponds to the contribution of 
the blue and green line segments xB and xF in Fig. [3] 
It receives no contribution from the "diagonal" region 
^,f(£i£) (the red line segment xD in Fig. [3j) , the corre- 
sponding contributions cancel out between the first and 
the second terms in Eq. (|4"3"|) (to the Q(l/N c ) accuracy). 
By this reason the conclusion in Ref. [27| that the evo- 
lution of the function T q .F{x, x) mainly follows a pattern 
determined by the quark splitting function can be mis- 
leading. 

The evolution does simplify, however, in the large x 
limit in which case the integration regions shrink to a 
point. One obtains 

n±T q , F {x,x) = ^ f d lp^-\ z )r q ALO , (47) 

dfi 7T J x £ 

where, retaining singular terms at z —> 1 only 



f 



(1-z). 



-6(1 -z) 



-N c S(l-z). 
(48) 



This result can be compared with the evolution of the 
usual Fi(x,Q 2 ) structure function which involves, to the 
same accuracy 



qq 



(z) = 2C f 



(49) 



and the twist-three contribution to the structure function 

9 2(x,Q 2 ) gam 



(1-z). 



(50) 

The last term in 5(1 — z) is written in the large-7V c 
limit. The contributions ~ 1/(1 — z) + are the same in 
all three cases, which indicates that all three functions 
T q , F (x,x; Q 2 ), Fi(x,Q 2 ) and g t 2 w '^(x, Q 2 ) may have the 
same functional dependence on the Bjorken variable x in 
the x — > 1 limit. Different terms ~ 8(1 — z) suggest, on 
the other hand, that twist-three functions are suppressed 
at large scales Q 2 compared to the twist-two distribution 
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1~ g ,F(x, x; Q 2 )/F 1 (x, Q 2 ) 
ff| u '- 3 (x,Q a )/F 1 (x s Q a ) 



a s (Mo) 

a s (Q) 
a s (Mo) 



2N c /b 



N c /b 



(51) 



The suppression of (x, Q 2 ) compared to F\{x,Q 2 ) 

exactly corresponds to the gap between the lowest 
anomalous dimension in the spectrum of twist-three oper- 
ators and the usual twist-two anomalous dimension. For 
the function T qt F(x, x; Q 2 ) we predict a stronger suppres- 
sion which translates to scaling violation in SSA. This 
result can be phcnomcnologically relevant. 



IV. FLAVOR-SINGLET EVOLUTION 

In the flavor-singlet sector one has to take into account 
mixing between the quark-antiquark-gluon and three- 
gluon operators with the same C-parity. Namely, 6 + 
gets mixed with T + and 6~ with T~ . For each case, 
the evolution equation takes the matrix form 



where 



l FQ H FF 



(52) 



(53) 



In what follows we assume that the quark-antiquark- 
gluon flavor-singlet operator is defined including the sum 
over rif light flavors, qFq — uFu + dFd + . . .. 
For the diagonal entries we obtain 



: + 4n f n 



13 I 



where H is given by Eq. (|37|) and 

B$ F = N c (n 12 + W 23 + Wsi - 4(W+ + n+ 3 ) 
- 2(H+ + H+) ± Q{Ui2 + Huj) - b , 



(55) 
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with &o = —N c — -nf. The off-diagonal entries in coor- 
dinate space take the form 



% F = -tn/zia Wa + «13 T 2Wr 3 . 



l FQ =iN c (l-P 23 ) — 

%12 



FQ 



N 2 - 4 
-t— (I + P23) 



2W+P 13 + 1 

1 r 



Zl3 



n , 

2W+ P13 - 1 



(56) 



The corresponding expressions in momentum space are 
given in the Appendix, Eqs. (|AT25l) - (fO0| . 

The evolution equations for conventional T q p q , AT q p q , 
T^p(x) AT^p can readily be obtained from Eq. ([B"2"|) by 
symmetrization in the arguments, as specified in Eq. (|3ip . 



In the limit of zero gluon momentum our result for the 
evolution of T qi p(x,x) = T q p q (— x, 0, x) differs from the 
corresponding expression in Eq. (107) of [13] by the same 
two terms as in the nonsinglet case; the terms in P qg (in 
our calculation due to Hq^) coincide. 

To compare our results for three-gluon distributions we 
introduce the functions 



Tp [x 1 x ) — ^ -^3f( x j x x , x) , 
ATp(x,x') = ~AT± F (-x',x-x',x), (57) 



which coincide with 



Tg,f( x , x') 



T, 



(d) 



(54) T { A f l F (x,x'), T^ F (x,x') defined 



(d) 



1 g,f( x , x ), 
in [27| , respec- 
tively. We remind that the AT F (x, x') distributions 
and hence T^q F (x,x'), T^q F (x,x') can be expressed 

in terms of T F {x,x') alias T G d F (x,x'), Tq' f {x,x'), so 
that they do not need to be considered separately. After 
some algebra one obtains the following equations for 
Tp(x,x): 



^T±(x,x) 



a s N r 



z 



-ii-:>(* + ~)r±(eo 



l + z 



1 - z 
± 



[T±«,x)-T±(€,0] 



T±(x,0-ATp(x,0 



T 3(1 - 2) [T±(x,x - - AT±(x, x - 



^p^) k f (£, ± t 9 , f (-£, -£) 



(58) 
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where 



and 



PM = 



A+ = l, 



(59) 



z 1 - z bo ... . 



l + (l-z) 2 



(60) 



Our result does not agree with that of Kang and Qiu 
[27j |. Eqs. (109), (HO) in an overall sign in front of the 
T AGF distribution (which may, however, be an artifact 
of a different sign convention for e^ ua f}), and two extra 
terms (up to the a s N c /i: factor) 



T±(x,x)± f 

J X 



(T*(x,x-Z)-AT*(x,x-£)) 

(61) 



which seem to have the same origin as the extra contribu- 
tions that we also have for the flavor-nonsinglet distribu- 
tions: The first term in (|61[) originates from contributions 
of the type in Eq. I]45p that involve a subtlety in taking 
the X2 — » limit, and the second one corresponds to the 
contribution of "exchange" type diagrams, the kernel Ti-i 2 
in (|55p , that give rise to mixing of regions with different 
partonic interpretation. 

Closing this section, we want to stress that the scale 
dependence has to be studied using complete evolu- 
tion equations for the three-particle parton distributions, 
Eq. (|52|) . Using the gluon-pole projected equations f58j) 
with a certain ansatz for the "off-diagonal" correlation 

functions Tqp (x, x') = T^ F (—x', x' — x, x) may be mis- 
leading as they are modified by the evolution themselves. 
We note in passing that the particular ansatz proposed 
in HI 



T</f(x 1 ,x 2 ) 



1 



T^\x uXl ) 
-T</f(x 2 ,x 2 ) 



-(x 1 -x 2 ) 2 /2a 2 



(62) 



and T^(fp(xi, x 2 ) — AT^ F (— x', x' — x,x) — 0, is incon- 
sistent with the constraint (|3"2")) . 



V. CONCLUSIONS 

We have given a complete reanalysis of the scale de- 
pendence of twist-three three-particle correlation func- 
tions that are relevant for calculations of single transverse 
spin asymmetries in the framework of collinear factoriza- 
tion. The calculation is done using the two-particle ker- 
nels for the renormalization of light-ray operators in the 
spinor basis, which are available from Ref. [351 ] . Evolution 



equations are derived for arbitrary parton momentum 
fractions, for the flavor-nonsiglet quark-antiquark-gluon 
distribution, Eqs. ([35|l . (|37|) . ([38| . and for the mixing 
matrix of the flavor-singlet quark-antiquark-gluon and 
three-gluon distributions with both positive and negative 
C-parity, Eqs. f52] ) -([55 ]) and (jX2l^ - (|X30"|) . Specializing 
to the case of zero gluon momentum we have compared 
our results with the recent calculation in Ref. [13] • There 
are two terms where we disagree, and their origin could 
be identified. As a byproduct of our calculation we pre- 
dict logarithmic scaling violation in the SSA at large val- 
ues of Bjorken x which may be phcnomcnologically rel- 
evant. Numerical studies of the evolution effects on re- 
alistic models of parton distributions will be considered 
elsewhere. 
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APPENDICES 

APPENDIX A: BFLK KERNELS IN 
MOMENTUM REPRESENTATION 

Gauge-invariant A-particle quasipartonic light-ray op- 
erators can be defined as a product of "plus" fields 

0(zi, ...,z N ) = CX(zi) ® . . . ® X(z N ) ee 
= Ci,... iN ([0, z 1 ]X(z 1 )) il . . . ([0, z N }X(z N )Y N ,(A.l) 

where X(z k ) = {tp +) $ + ,x+,X+, /++, /++}, [0,2fc] are 
Wilson lines in the appropriate representation of the 
gauge group, ii, . . . , ijv are color indices and Ci 1 ...i N is 
an invariant color tensor such that 

[(*l)*i»i + (*2)fe 2 j 2 + • • • + (tN)k N i N ] S ii,-,iN = • 

(A.2) 

Here and below it is implied that the generators t a are 
taken in the appropriate representation, 

(t a ^) 1 = T^ip 1 ' for quarks ip, \ 
t a X = { {t a ^Y = -Tfaip*' for antiquarks \ 
(t a f) b = i i bab 'f for gluons /, / 

(A.3) 

where T a are the generators in fundamental repre- 
sentation. The condition in Eq. (|A.2[) ensures that 
0(zi, . . . , zm) is a color singlet. 
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For each N, the set of quasipartonic operators with the 
same quantum numbers is closed under renormalization 
[3ll | . A renormalized quasipartonic operator is written as 



[Oi(X)] R = Z tk O k (X ) 



(A.4) 



where X = ZxX is the bare field. Renormalized opera- 
tors satisfy the RG equation 



; ••■ \ (A..-. 1 



Here (3(g) is the (QCD) beta function and 



7 = -u— ZZ _1 



is the matrix of anomalous dimensions. To the one-loop 
accuracy one obtains in dimensional regularization (D = 
4- 2e) 



Z = I- 



47T£ 



and 



2tt 



(A.6) 



The operator H (Hamiltonian) is given by the sum of 
two-particles kernels, H- fc 



N 
i,k 



(A.7) 



The general structure of the kernels is 

m^ 2) [x^( Zl )®x l -(z 2 )] = 

= EEi c ^^ ri ®^(^^)- (A- 



q i[i' 2 



Here [£?,,]£*? is a color tensor, ?4f is an SX(2,K) invari- 
ant operator which acts on coordinates of the fields, and 
q enumerates different structures. 

Explicit expressions for the two-particle (BFLK) ker- 
nels in the light-cone gauge in coordinate space are given 
in Table U [35| . We tacitly assume existence of rif quark 
flavors; the kernels and B s correspond to the flavor- 
nonsinglet and flavor-singlet quark pairs, respectively. In 
the kernel B s the generators t a should be taken t a = T a 
for the ip ® ip pair, and t a — (— T a )' for the x ® X pair. 

The kernels in Table U are written in terms of several 



"standard" SX(2)-invariant operators defined as 

[H (p](z 1 ,z 2 ) = [ — 2ip(zi,z 2 ) 
Jo " L 

-a 2 ^-\(z^ 2l z 2 )-a 2 ^-\(z u zq i ) 
-1 



[H d <p]{Z!,Z2) 
[H + <p]{Z!,Z2) 

[H+v](zi,z 2 ) 







dttfi*'- 1 ^- 1 ^,^), 

'da[ a dpa 2 ^- 2 ^- 2 ^ 2l z^), 



a(3 



Jo 

1 />a 

da / d/3a 2jl ~ 2 /3 2j2 " 2 . 
Jo \<xp 

x ^"2,^21) > 



[7i-^(2i,2 2 ) = / da d(3a 231 ' 2 (3 2l2 - 2 ip(z^ 2 ,z^). 

Jo J a 

[H e 12 { %}( Zl ,z 2 ) = [ daa 2h - k - 1 a k - 1 cp(z? 2 ,z 2 ), 



(A.9) 

where in the last line it is assumed that < k < 2j\. 
Here and below we use a shorthand notation zf k = Zi6t + 
z k a , a = 1 - a, Z\ 2 = z\ — z 2 . 

The kernels depend on the conformal spins of the par- 
tons they are acting on, j = 1 for quarks and antiquarks 
and j = 3/2 for gluons. Another shorthand notation is 

J a (zi,z 2 ) = 
= ]T (yi(zx)T^i(z 2 ) + xi( Zl )T a x i(z 2 )) , 

A 

J ac (zuz 2 ) = 
= ]T [i>i( Zl )T a T^i(z 2 ] - xi(z2)T c T a xi(zi) 

A 

(A.10) 

where A is the flavor index; the sum runs over all possible 
flavors. The constants 



b /4,N c , 



(A.ll) 



3 

a q = - , a 

where 60 = if^c — correspond to the "plus" quark 
field and transverse gluon field renormalization in the 
axial gauge 



Z q — 1 + - — o„Cf 
Zire 



Z g = l + ^o g C A . (A.12) 

Z7T£ 



Finally, P\ 2 and P ac stand for permutation operators in 
position and color space, respectively. For example 



PxiJ ac (z Xl z 2 ) = J ac (z 2 , Zl ). 

PacJ aC (zi,Z 2 ) = J Ca ( Zl ,Z 2 ). 



(A.13) 



Going over from the light-ray operator renormaliza- 
tion to evolution equations for parton distributions cor- 
responds to the Fourier transformation of the functions 
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i 




x 2 ] 


A 

A 


V>+®i/>+, V>+®V>+, 
V>+®x+, x+®x+, x+®x+ 


-2(t a .,t a ;)\n- 2a 


,]x l i(zi)®X l 2( Z2 ) 


B NS 


^+ ® x+,Vm- ® x+, 
ip+ ®V>+,x+ ® x+ 


-2(^ 1!; ^)[H-W + -2a, 


® X*=(z 2 ) 


B s 


Vm- ® x+ ®X+ 


-2(*? lii *? a i)[«-W + -2ff," 
-41% H d J a {z 1 ,z 2 ) 
-2iz 12 {{t a t b ) tJ \n + +H+] 4 


2(^ a ) lJ H-}/^ + ( 2 i)®/i + ( 22 ) 


c 




-2(C'®4') [W-<7,-<7 9 ] 

-2(t°'*°) ii /ft a «Si (1) X" , (« 1 ) 


X a \ Zl ) ® X 1 ' (z 2 ) 
®X l \z 2 ) 


D 


n + ®^+, n+®x+ 


-2(4,,® 4,) [£-2H+-a 


q - a a\ 


X a '( 2l )®X l '( 22 ) 


E 




~~ 2(^aa'*cc0 


W-2a 9 


X a ' ( Zl ) ® X c ' (z 2 ) 


F 


n+ « /++ 


— 2(t aa 't C c') 

+i2(iL'4': 

4— [2W+P 

Z12 | 


12 - Pac] 


- 2 H + - 2a g 

n) ® f++(za) 
\ -6"ft d ) J ac 


]/?'+(») 

(Z2,2l) 



TABLE I: Summary of the BFLK kernels 



n d 7i e,(k)^ n+ n+ n - from coor dinate to the mo- 
mentum fraction representation. 

A generic two-particle kernel in momentum space has 
the form 



[Hep] (ui,u 2 ) 



T>vH(u 1 ,u 2 \v 1 ,v 2 ) tp(vi,v 2 ) , 

(A.14) 



where T>v = dvidv 2 S(ui+u 2 — vi — v 2 ). In this Appendix 
we use the letters u, v for momentum fractions. Note that 
the integration regions over v\ and v 2 in Eq. (|A.14[) are 
(formally) infinite; in practice they are constrained by 
the support properties of parton distributions and also 
(Heaviside) step-functions that are present in the kernels. 
For later convenience we introduce a notation 



e(oi, . . . , on) = n AM - n ^(- a <) ■ ( a - is ) 



k=l 



k=l 



We remind that the kernels depend on the conformal 
spins of fields they acts on, cf. Eq. (|A.9|) . so that, in 
principle, each one has to carry a pair of indices (ji,j 2 ) 
which we omitted for brevity. 
We obtain 



nUu\v) = 



Q(ui,u 2 ) ( u 



Ui+U 2 \Ui+U 2 



«2 



Ui + U 2 



(A.16) 



Htf\u\v) = 6( " 1 '" 1 "" l) ( ^-) 2n ~ k ( 1 - ^iV" 1 , 



(A.17) 



2ii-2 



= ^KgMt-ijO y( tt )-fa) <p{v) 



u 2 9(w 2 , v 2 - u 2 ) 
v 2 v 2 - u 2 



u 2 \ 



2J2-2 



r(«)-(^J *>(«) I > (A.18) 



Note that Hf 2 (|A~16|) vanishes in the DGLAP region. 

For the remaining kernels the expressions for arbitrary 
conformal spins become too lengthy, so we specialize to 
the particular cases of interest. 

The operator Ti + enters the BFLK kernels in Ta- 
ble U with the values of conformal spins (ji,j 2 ) = 
{(1, 1), (3/2, 1), (1, 3/2), (3/2, 3/2)}. In momentum frac- 
tion space one obtains 



H + ' {h ' h) {u\v) 



>h)i 



- Q(—Ui, u 2 , u\ — v^A^ 1 
+ e(ui,ua,«a - u 2 )B^' j ^(u\v) 

(A.19) 



The first term contributes to the DGLAP region only, 
the other two — to the ERBL region. 
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One derives the following expressions: 



,4(M) 
5(1,1) 

Ad' 1 ) 
Bd.i) 



1 1 

- 2 ' 2 



'3 3 
V 2 ' 2 



u 



u 



u 



vi - Ui 

ViV 2 
U 2 



v 2 (u 1 + u 2 ) 

vi(ui + u 2 ) ' 
1 v\ — u\ 



2 v\v 2 ' 

1 u 2 u 2 + 2ui 

2 v 2 (ui + u 2 ) 2 ' 
luf v 2 + 2v x 
2^( Ul +u 2 y 

1 Vi — Ui 



2 v\v\ 



1 fU2 

2 \v 2 J ui +u 2 

2 



WlU 2 + UiV 2 - -(«! - Ui) 2 



1 



1 /ui 

2 Ui 



1 



Ul + u 2 



UiV 2 - U 2 Vi/?> 

(mi + u 2 ) 2 

u 2 vi — uiv 2 /3 
(ui + u 2 ) 2 

(A.20) 

Obviously 

n + ^ 2 Hui 1 u 2 \vi,v 2 ) = n+^ 2 - 1 Hu2,ui\v 2 ,vi). 

The "modified" Tt + kernel is needed for the conformal 
spins (ji,j 2 ) = (3/2,3/2) only. It is given by a similar 
expression 

H+>Vi>M(u\v) = e{-ui,u 2 ,ui - vi)A {n ' l2 \u\v) 

+ e(ui,u 2 ,v 2 - u^B^ 1 '^ (u\v) 

+ Q(u 1 ,u 2 ,v 1 -ui)& n ' h) {u\v) 

(A.21) 



with 



A(u\v) 
B(u\v) 
C{u\v) 



1 (Ui - Vi) 3 

6 v\v\ 

1 fu 2 \ 2 UlV 2 - Ulti 2 /3 

2 \v~2~J {ui+u 2 f ' 

1 / Ui\ u 2 vi - v 2 Ui/3 

2 \v^J ( Ul +M 2 ) 3 



(A.22) 



Note that the factors A, B, C (and, similar, A, B, C) sat- 
isfy the relation 

The kernels 7i + (it|i>) and 7Y + (tt|t;) are continuous func- 
tions in the whole domain although they are given by 
different expressions in the DGLAP and ERBL regions. 

The kernel Tt~ has a different structure of regions, 
namely 

H~(u\v) — 9(ui, —u 2 ,u 2 — vi)D(u\v) 
+ <d(ui,u 2 , v 2 - ui)E(u\v) 
+ Q{ui,u 2 ,ui-v 2 )F(u\v). (A.23) 



For the spins of interest (ji,j 2 ) = {(3/2, 1), (3/2, 3/2)} 
one gets 

2 



F (|,l) 

F (|,|) 



1 / ^2 - «1 

2v 2 \ w i 
J_ / ui 
2v 2 \ui+u 2/ 
u 2 2u\Vi — u 2 v 2 
2v} (ui+u 2 ) 2 ' 

1 (U 2 - Vl) 3 

6 v\v\ 

1 /mi\ 2 U 2 V 2 — Ullil/3 



2 \v 2 J (ui + u 2 ) 3 

1 ( u ^\ 2 u i w i — v 2 u 2 /i 

2 V^J (lii+M 2 ) 3 



(A.24) 



We also give here the explicit expressions for the off- 
diagonal kernels (|56[) in momentum representation. The 
Hg^-kernel can be written as 

M± F (u\v) = n f (V? 3 (u\v) T Vi 3 (u\v)) , (A.25) 

where V^~ 3 (u\v) and Vi 3 (u\v) have the decomposition 
over different regions of the form (|A.19|) and (|A.23[) , re- 
spectively. For the corresponding functions Ay,i?y,CV 
and D\>, Ey, F\> one obtains 



A v (u 
B v (u 
C v (u 
D v (u 
E v (u 
F v {u 



UlU 3 
2 2" ' 

UiU 3 Vi + 3v 3 
v 3 (ui + u 3 ) 3 ' 

uiu 3 v 3 + 3vi 
vj (ui + u 3 ) 3 ' 

{u 3 - Vl) 2 

2 2 ' 



v 3 (ui + u 2 ) 2 
v{{ui +u 2 ) 2 



2u 3 v 3 



ui{ui + u 3 ) 
2uiVi 



u 3 (ui + u 3 ) 



(A.26) 



The H| Q -kernel can be written as 



I+ Q = JV C (1 - P 23 ) [W+ +W-- 2AW] , 

A 2 - 4 

IT - - 



-(I+P23) [w + + w- 



(A.27) 



Here 



AW(w|u) = 8(113, "3 ~ u 3 ) - Q(ui,u 3 ) 



u\{Zu 2 + Ml) 

(ui + u 3 ) 3 

(A.28) 



The kernel W + (ui , u 3 \vi , v 3 ) has the structure (|A.19|) 
with 

A w (u\v) = l, B w {u\v) = -C w {u\v) = ^, (A.29) 
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and the kernel W~ (ui,us\vi, V3) takes the form (|A.23|) The kernels W ± ,A>V are antisymmetric under 
with (ui,Vi) <-> (u2,V2). Note that the kernels in (|A.27|) have 

discontinuities on the DGLAP-ERBL boundaries. 

D w (u\v) = , 



E w (u\v) 



1 



2 v 3 (ui + u 2 ) ' 



1 u 



2 



F w {u\v) = -^- + — — 2 (A.30) 

v 1 {u 1 +u 2 ) 
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